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Abstract 

^ ■ In this paper we employ the hydrodynamic formulation of time-dependent density functional 

^ theory to obtain the van der Waal coefficients Cq and Cs of alkali-metal clusters of various sizes 

<Z3 ■ 

.^including very large clusters. Such calculation becomes computationally very demanding in the 

. 

^orbital-based Kohn-Sham formalism, but quite simple in the hydrodynamic approach. We show 

i ^ that for interactions between the clusters of same sizes, Cq and Cg sale as the sixth and the eighth 

t— I power of the cluster radius rsepectively, and approach the respective classically predicted values 

r~- for the large size clusters. 
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I. INTRODUCTION 



The long-range van der Waals forces play an important role in the description of many 
physical and chemical phenomena such as adhesion, surface tension, physical adsorption, etc. 
The correlations between electron density fluctuations at widely separated locations give 
rise to these long range forces. For clusters the knowledge of coefficients of van der Waals 
interaction is useful in describing the cluster-cluster collisions and also for characterizing the 
orientation of clusters in bulk matter. The van der Waals interaction coefficients, namely 
Cq B and C% B between small alkali metal particles within the spherical jellium background 
(S JBM) model of have been theoretically calculated in the past using time dependent 

density functional theory (TDDFT) [4] . In Refs. time-dependent Kohn-Sham (TDKS) 

4] formalism of TDDFT was employed to obtain C$ B and C$ B . On the other hand, in Ref. 

3|], we employed density based method within TDDFT formalism to obtain the van der 
Waals coefficient Cq B for the clusters. The TDKS formalism is an orbital based theory and 
involves calculation of these orbitals in a self-consistent manner. Thus, as the size of system 
increases these calculations tend to become computationally cumbersome. In such cases, 
our approach [3J] is much easier to as it makes calculations size independent with hardly any 
loss in accuracy. 

The mam aim of this paper is to extend our previous study fl El by applying the density 
based hydrodynamic approach to calculate the higher-order van der Waals coefficient C$ B 
of alkali atom clusters characterizing dipole-quadrupole interaction of long-range force and 
study its evolution with the size of clusters. In addition to this we also calculate C$ B and 
C^ B coefficients for the pair interaction between potassium clusters of various sizes and 
between potassium and sodium clusters of different sizes as well. 

Before proceeding further, it is necessary to note that the density functional theory (DFT) 
in principle should give the exact ground-state properties including long range van der Waals 
energies. However, the frequently used local density approximation (LDA) [3, 0] and gen- 
eralized gradient approximations (GGA) (3, [lO, [ljj fail to reproduce the van der Waals 
energies. This is due to the fact that LDA and GGA cannot completely simulate the corre- 
lated motion of electrons arising from coulomb interaction between distant non overlapping 
electronic systems. It is only recently that several attempts 13, 3] have been made to 
obtain van der Waals energies directly from the ground-state energy functional by correcting 
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the long range nature of the effective Kohn-Sham potential. On the other hand, it is possi- 
ble to make reliable estimates of the van der Waals coefficients by using expressions which 
relate these coefficients to the frequency dependent multipole polarizabilities at imaginary 
frequencies 15, 3]. We follow this route for the calculation of these coefficients. 

The paper is organized as follows: In section II we express Cq and C§ in terms the 
dipole and the quadrupole dynamic polarizabilities. We then briefly describe our method 
of obtaining these polarizabilities employing hydrodynamic approach of TDDFT. Results of 
our calculations are presented and discussed in section III. 

II. METHODS OF CALCULATION 

The nonretarded electrostatic interaction energy between two spherically symmetric elec- 
tronic systems A and B separaterd by interaction distance R can be written as [15I . 1|| 

riAB riAB 

The above expression has been obtained by assuming that R is very large so that the charge 
distributions corresponding to the two systems do not overlap. The coefficient Cq B describes 
the dipole-dipole interaction, and C$ B corresponds to the dipole-quadrupole interaction 
between system A and B. These coefficients can be expressed in terms of the dynamic 
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arizability ai(u) (where I denotes the index of multipolarity) by following 



18|: 



C£ B = C(A,1;B,1) 

C£ B = C(A,1;B,2) + C(A,2;B,1) (2) 



C(A,h;B,l 2 ) = } 2 !' + ^\. I duai{iu)a B {iu) (3) 



with 

r,( A L ■ R 1„\ = 

2tt(2/ 1 )!(2/ 2 )! J 

where af(iuj) is the multipole polarizability of system X (X = A or B) at imaginary 
frequency u = iu. Although af(iu) does not have any physical significance, expressions 
given by Eqs. (J2J) and (JHJ) make the computation of C§ B and C% B straightforward. Moreover, 
mathematically af(iu) is better behaved than its real frequency counterpart a(u): it does 
not have any singularity and decreases monotonically from its static value af(0) to zero as 
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uj — > oo. Consequently, the quadrature in Eq. Q can be computed quite accurately. To 
determine the frequency dependent polarizabilities appearing in Eq.©, we employ variation- 
perturbation method within the hydrodynamic approach of TDDFT. We now describe the 
theory in brief. For details the reader is referred to the literature 

The basic dynamical variables of the hydrodynamic theory are the time dependent density 
p(r, t) and the velocity potential S(r, t). Thus the total time-averaged energy can be ex- 
pressed in terms of these two variables. For our purpose we need to evaluate the second-order 
change in the time-averaged energy as this is directly related to the frequency dependent 
multipole polarizability by the relation 

ai (u) = -4£f (4) 

(2) 

The second-order time-averaged energy E\ } in turn can be expressed as 

+ J gg^V >(r,t)rfr + I J (VS"> ■ VS"V>(r)<fcj , (5) 

where the curly bracket denotes the time averaging over a period of the applied oscillating 
field and p(°)(r) represents the ground-state density. It is easily shown [3] that £j is 
stationary with respect to the variations in the first-order induced density p^(v,t) and 
the induced current-density S^(r,t). Consequently, £j can be determined by choosing 
appropriate variational forms for p^\r,t) and S^(r,t) and making E^ stationary with 
respect to the variations in the parameters of p^(r, t) and S^(r, t). In the above expression 
the functional 

E[p}=T s [p}+E H [p}+E xc [p}, (6) 

where T s [p], E H [p] and E X c[p] denote the kinetic, Hartree and the exchange-correlation 
(XC) energy functionals respectively. The exact forms of T s [p] and Exclfi] are n °t known. 
Consequently to perform any calculation one needs to use approximate forms for these 
functionals. On the other hand, the hartree energy functional En[p] representing classical 
coulomb energy is exactly known and it given by 

EM = \JJ (7) 

For the purpose of calculation the multipolar applied potential v l app (r,t) is chosen to be 

v l app (T,t) = Sr l Y l0 (e,<f>)cosut (8) 



where 8 and u> represent the amplitude and the frequency of the applied periodic electro- 
magnetic field. In accordance with the above form of the applied potential the variational 
forms for p^'(r, t) and S^'(r, t) are chosen to be 

p^fr, t) = p^\r , u)cosujt 

S {1 \r,t) = uS {1) {r,u)sinujt (9) 

with 

i 

S (1) (r,w) = J>rXr)^ (M) (10) 

i 

where p^ (r) is the ground-state density and C; and di are the variational parameters obtained 
by minimizing time-averaged second-order energy . On substituting Eq. © in Eq. (jSJ) 
and taking average over time we get 

+ ^ sW(r,u)pW{r,u;)dr+ ( £ J (VS« • VS^)p^(r)dr, (11) 

At this point it is important to point out that the VP method discussed above is also 
applicable for the imaginary frequencies with uj 2 replaced by — u 2 in Eq. (fTT|) [3] . This allows 
us to determine dynamic multipolarizability at imaginary frequencies (a(iu)) by exactly the 
same procedure as employed for getting a(u). All that is required for this is to change uo 2 
to — uj 2 in Eq. (Jllj) . This is done very easily in the numerical code written for determining 
dynamic polarizability at real frequencies. 

As mentioned earlier the calculation of E\ requires app roximating the functionals T s [p] 
and Exc[p\- To this end we choose the von Weizsacker [l9j form for T s [p] which is given as 

Twlp] = lJ^l d , (12) 

Our previous experience with the calculation of response properties has led us to choose von 
Weizsacker KE functional for the polarizability calculation. For the XC energy, adiabatic 
local-density approximation (ALDA) [4] is accurate enough to describe the energy changes. 



5 



Thus the exchange energy is approximated by the Dirac exchange functional 

E x [p] = C X J pl(r)dr 



C x = --A~) • d3) 



4 \H ^ 

and for the correlation energy functional we employ the Gunnarsson-Lundqvist (GL) |2l| 
parametrized form for it . 

In the present paper the ground-state densities p^(r) of clusters are obtained by em- 

nn 

ploying purely density-based extended Thomas-Fermi (ETF) [22J, |23j method within the 
spherical jellium background model (SJBM) of metal clusters. This approach yields the 
ground-state densities of very large clusters (containing up to 1000 atoms) easily thereby 
allowing us to study the evolution of van der Waals coefficients with the size of clusters. 
For details of the ETF method and its application to study alkali-metal cluster, we refer the 



reader to 



23, 



24j. In the next section we present the results for Cq B and C$ B of 



alkali metal clusters by employing the method describe above. 



III. RESULTS AND DISCUSSION 

We have performed calculations for the coefficients C§ B and C$ B between clusters of 
alkali metal atoms of various sizes. In this paper we consider clusters made up of sodium 
(r s = 4.0a. u.) and potassium (r s = 4. 86a.it.) atoms ( where r s is the Wigner-Seitz radius of 
the cluster species). First we discuss the results for the sodium clusters. For completeness 
we also include Cq B numbers of sodium clusters in this paper. As in our previous study, 
Table I and II we present the results for C$ B and C$ B respectively, between sodium clusters 
containing 2, 8,20, and 40 atoms. For comparison we also show the corresponding TDKS 
results of Ref. |2j in parenthesis. Table I and II clearly show that the numbers obtained by 
employing the hydrodynamic approach for small clusters are quite close to the corresponding 
numbers obtained by the TDKS method. The numbers for Cq B obtained in this paper are 
systematically lower than the orbital-based TDKS results. The difference between the two 
results is slightly more for larger clusters than the smaller ones, the maximum difference 
being of the order of 10%. Now we discuss the main results of this paper, that is, numbers 
for Cg B between clusters of sodium clusters. 
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In Table II we present the results for the coefficient C§ B . This table clearly shows that our 
numbers for C% B are quite close to the corresponding TDKS results. In comparison to the 
results for Cq B our numbers for C$ B are closer to the corresponding TDKS values. Moreover, 
unlike Cq B case, the numbers for C$ B are not always lower than the corresponding TDKS 
results. The values of C$ B between 2 and 2 and, 2 and 8 atom sodium clusters obtained by 
the hydrodynamic method are slightly more than the corresponding TDKS results. We next 
apply hydrodynamic approach to calculate the van der Waals coefficients between potassium 
clusters of various sizes. 

In Table III and IV we present the results for Cq B and C$ B respectively, between potas- 
sium clusters of different sizes. Again two tables clearly show that both Cq B and C$ B 
between potassium clusters are obtained quite accurately with the hydrodynamic approach. 
The numbers obtained by us for potassium clusters are all lower than the corresponding 
TDKS results. Table III and IV also clearly show that the difference between the hydro- 
dynamic and the TDKS results is lower for the larger clusters than than the smaller ones. 
On the other hand, for sodium clusters the difference between the two results shows the 
opposite trend. 

Next we now present in Table V and VI the results for Cg and C AB respectively, for 
the pair interaction between sodium and potassium clusters with 2, 8, and, 20 constituent 
atoms. Once again we see from these two table that hydrodynamic numbers quite close to 
the corresponding TDKS results barring few exceptions like the number for Cq B between 
two types of clusters containing 8 atoms each. Table V and VI also show that like the pair 
interaction between clusters of identical atoms the value of C$ B for two different clusters 
obtained by hydrodynamic approach is closer to the TDKS as compared to the corresponding 
Cq B results. 

With favourable comparison between the numbers obtained by our approach and those 
by the TDKS approach for small clusters we now employ the hydrodynamic approach to 
evaluate the van der Waal coefficients for much larger clusters. This allows us to study the 
size dependence of van der Waals coefficients and how they approach the asymptotic classical 
values. In the present paper we have performed calculations for clusters containing up to 
1000 atoms without any increase in computational effort. To study the size dependence of 
van der Waals coefficients we plot in Figs. 1 and 2, Cq B / 'R® and C$ B / 'R* respectively, as 
functions of R ( where R = TgN 1 / 3 and N denotes number of atoms). These two figures 
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clearly exhibit that as the size of cluster increases the values of coefficients C aa /Rq and 
C$ a /Rq saturate to constants numbers, indicating that C AA and C AA scale as the sixth and 
the eight power of the radius Rq of the cluster respectively. This trend is consistent with the 
fact that the properties of metal clusters approach their corresponding classical values as 
the size of cluster is increased. The classical expressions for C AA /R® and C AA /R$ between 
clusters of same sizes can be written as 

Cq A /Rq = -u Mie 

C£ A /R 8 = H—^—^. (14) 
Here uJuie is the Mie resonance frequency given by 

&Mie — \ ~ (15) 
y 's 

and it is equal to 1/V3 bulk plasmon frequency. The classical expressions for C aa /Rq and 
C$ a /Rq given above are derived in Ref. |2j by assuming that all the strength of respective 
multipole resonance is concentrated in a single peak. By substituting the values of r s in 
Eq. (|I1J) we get following numbers for sodium clusters C AA /R^ = 0.094 and C AA / R% = 0.49 
whereas for potassium clusters Cq A /R® = 0.07 and C AA /Rq = 0.37. These classical numbers 
for C AA /Rq and C AA /R$ are shown by asymptotic straight lines in Figs. 1 and 2. We see 
from Figs. 1 and 2 that the hydrodynamic approach yields correct asymptotic values of the 
van der Waals coefficients or the values to which these coefficient saturate as the size of 
cluster grows. 

To conclude, we have extended the applicability of the hydrodynamic approach within 
TDDFT to the calculation of higher order van der Waals coefficient C$ B between clusters 
of various sizes and different species. Our results for both C AB and C AB are quite close to 
more accurate orbital based TDKS approach. In particular we have found that the numbers 
obtained by hydrodynamic approach for C AB are more accurate than that of C AB . For both 
sodium and potassium clusters we have been able to calculate C$ B and C AB coefficients for 
clusters containing up to 1000 atoms. Thus we have been able to get the evolution of these 
coefficients as a function of cluster size and have shown that they approach their respective 
classical values for large clusters. 
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TABLE I: Dispersion coefficient C% for sodium atom clusters in atomic units (a.u.). The numbers 



follow the notation 3.60(3) = 3.60 x 10 3 . The numbers in parenthesis are results of Refs. 



N 


2 


8 


20 


40 


2 


2.56(3) 


9.47(3) 


2.27(4) 


4.42(4) 




(2.62(3)) 


(1.02(4)) 


(2.45(4)) 


(4.74 (4)) 


8 




3.51(4) 


8.42(4) 


1.64 (5) 






(4.01(4)) 


(9.55(4)) 


( 1-86 (5)) 


20 






2.02(5) 


3.93 (5)) 








(2.29(5)) 


( 4.45 (5)) 


40 








7.63 (5) 










(8.60 (5) 



TABLE II: Dispersion coefficient Cg B for sodium atom clusters in atomic units (a.u.). The numbers 
follow the notation 3.60(3) = 3.60 x 10 3 . The numbers in parenthesis are results of Refs. 



N 


2 


8 


20 


40 


2 


3.08 (5) 


2.28 (6) 


8.95 (6) 


2.58 (7) 




(2.97(5)) 


(2.27(6)) 


(9.08(6)) 


(2.69 (7)) 


8 




1.26(7) 


4.33 (7) 


1.15 (8) 






(1.32(7)) 


(4.59 (7)) 


(1.25 (8)) 


20 






1.35 (8) 


3.37 (8) 








(1.44 (8)) 


(3.67 (8)) 


40 








8.00 (8) 










(8.82 (8)) 
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TABLE III: Dispersion coefficient C^ B for potassium atom clusters in atomic units (a.u.). The 

Q 

numbers follow the notation 3.60(3) = 3.60 x 10 3 . The numbers in parenthesis are results of Ref. |2| 

N 2 8 20 



2 5.35 (3) 2.03 (4) 4.93 (4) 

(6.28 (3)) (2.34(4)) (5.48(4)) 

8 7.70 (4) 1.87 (5) 

(8.71 (4) (2.02 (5)) 

20 4.55 (5) 



(4.74 (5)) 



TABLE IV: Dispersion coefficient C$ B for potassium atom clusters in atomic units (a.u.). The 
numbers follow the notation 3.60(3) = 3.60 x 10 3 . The numbers in parenthesis are results of Ref. 

B 



N 2 8 20 

2 8.69 (5) 6.87 (6) 2.77 (7) 

(1.10 (6)) (7.81 (6)) (3.03 (7)) 

8 3.97 (7) 1.38 (8) 

(4.37 (7)) (1.43 (8)) 

20 4.36 (8) 

(4.55 (8)) 
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TABLE V: Dispersion coefficient C§ B for the pair interaction between sodium clusters (values of 
N along the column) and potassium clusters (values of N along row) in atomic units (a.u.). The 
numbers follow the notation 3.60(3) = 3.60 x 10 3 . The numbers in parenthesis are results of Ref. 3] 



N 


2 


8 


20 


2 


3.68 (3) 


1.40 (4) 


3.40 (4) 




(4.01 (3)) 


(1.50 (4)) 


(3.49 (4)) 


8 


1.35 (4) 


5.16 (4) 


1.26 (5) 




(2.16 (4)) 


(8.05 (4)) 


(1.88 (5)) 


20 


3.25 (4) 


1.24 (5) 


3.00 (5) 




(3.75 (4)) 


(1.40 (5)) 


(3.26 (5)) 



TABLE VI: Dispersion coefficient Cg for the pair interaction between sodium clusters (values of 
N along the column) and potassium clusters (values of N along row) in atomic units (a.u.). The 
numbers follow the notation 3.60(3) = 3.60 x 10 3 . The numbers in parenthesis are results of Ref. 

a 



N 


2 


8 




20 


2 


5.21 (5) 


4.44 


(6) 


1.84 (7) 




(5.77 (5)) 


(4.41 


(6)) 


(1.84 (7)) 


8 


3.55 (6) 


2.23 


(7) 


8.29 (7) 




(4.52 (6)) 


(2.91 


(7)) 


(1.12 (8)) 


20 


1.35 (7) 


7.30 


(7) 


2.45 (8) 




(1.51 (7)) 


(7.75 


(7)) 


(2.57 (8)) 
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Figure captions 

Fig.lPlot of van der Waals coefficient C AA in units of of alkali-metal clusters: sodium 
(solid circles) and potassium (solid squares) as a function of Rq. The lines are drawn as a 
guide to eye and horizontal lines represent corresponding classical values of van der Waals 
coefficient. 

Fig.2Plot of van der Waals coefficient C AA in units of Rq of alkali-metal clusters: sodium 
(solid circles) and potassium (solid squares) as a function of i?o- The lines are drawn as a 
guide to eye and horizontal lines represent corresponding classical values of van der Waals 
coefficient. 
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